arXiv: 1504.07819v4 [math.PR] 4 Apr 2016 


EXTREMES OF THE SUPERCRITICAL GAUSSIAN FREE FIELD 


ALBERTO CHIARINI, ALESSANDRA CIPRIANI, AND RAJAT SUBHRA HAZRA 


Abstract. We show that the rescaled maximum of the discrete Gaussian Eree Eield (DGEE) 
in dimension larger or equal to 3 is in the maximal domain of affracfion of fhe Gumbel dis- 
fribufion. The resulf holds bofh for fhe infinife-volume field as well as fhe field with zero 
boundary conditions. We show that these results follow from an interesting application of 
fhe Sfein-Chen mefhod from Arrafia ef al. (1989). 


1. Introduction 

In this article we consider the problem of determining the scaling limit of the maximum 
of the discrete Gaussian free field (DGFF) on d > 3. Recently the maximum of the 
DGFF in the critical dimension d = 2 was determined in Bramson et al. (2013). In this case, 
due to the presence of the logarithmic behavior of covariances, the problem is coimected 
to extremes of various other models, for example the Branching Brownian motion and the 
Branching random walk. In d > 3, the presence of covariances decaying polynomially 
changes the setting but the behavior of maxima is still hard to determine (Chatterjee, 
2014, Section 9.6). This dependence also becomes a hurdle in various properties of level 
set percolation of the DGFF which were exhibited in a series of interesting works (Drewitz 
and Rodriguez (2015), Rodriguez and Sznitman (2013), Sznitman (2012)). The behavior 
of local extremes in the critical dimension has also been unfolded recently in the papers 
Biskup and Louidor (2013, 2014). 

We consider the lattice Z'^, d > 3 and take the infinite-volume Gaussian free field 
with law P on . The covariance structure of the field is given by the Green's 
function g of the standard random walk, namely E [(pa(p^] = g{a. — /3), for a, /3 G Z^. For 
more details of the model we refer to Section 2. It is well- known (see for instance Lawler 
(1991)) that for a; 7^ /3, g(a. — /3) behaves likes ||a; — and hence for ||a; — /3|| —^ -|-oo, 

the covariance goes to zero. However this is not enough to conclude that the scaling is the 
same of an independent ensemble. To give an example where this is not the case, when 
Vj^ is the box of volume N, XlaeW, order unlike the i. i. d. setting (see 

for example Funaki (2005, Section 3.4)). 

The expected maxima over a box of volume N behaves like ^yigiOjlogN. An indepen¬ 
dent proof of this fact is provided in Proposition 4 below; this confirms the idea that the 
extremes of the field resemble that of independent A/'(0,g(0)) random variables. In this 
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article we show that the similarity is even deeper, since the fluctuations of the maximum 
after recentering and scaling converge to a Gumbel distribution. Note that in d = 2 the 
limit is also Gumbel, but with a random shift (see Bramson et al. (2013, Theorem 2.5), 
Biskup and Louidor (2013)). The main results of this article is the following. 

Theorem 1. Let A he a subset with \A\ = N^. We define two sequences as follows: 


bN = 


so that for a// z G R 


/^^_loglogN + log(4;r) 


lim p <P,-bN 


N—>+00 y 


On 


and a]~^ = g{0){bfq) 


-1 


( 1 ) 


= exp(— e 


— Z> 


Note that scaling and centering are exactly the same as in the i. i. d. A/'(0,g(0)) case, see 
for example Hall (1982). As in d = 2, the argument depends on a comparison lemma. We 
show that in fact the proof is an interesting application of a Stein-Ghen approximation 
result. Not only does the result depend on the correlation decay, but also crucially on 
the Markov property of the Gaussian free field. We use Theorem 1 of the paper by Ar- 
ratia et al. (1989) which approximates an appropriate dependent Binomial process with a 
Poisson process, and gives some calculable error terms. In general showing that the error 
terms go to zero is a non-trivial task. In the DGFF case, thanks to estimates on the Green's 
function and the Markov property, the error terms are negligible. 

The techniques used for the infinite-volume DGFF allows us to draw conclusions also for 
the field with boundary conditions. For n > 0 let N := n‘^; we consider the discrete hy¬ 
percube Vjsf := [0, n — 1]^ n Z^. We define therein a mean zero Gaussian field (i/7fl;)^g^d 
whose covariance matrix i^))a,^€Viv is the Green's function of the discrete Lapla- 

cian with Dirichlet boundary conditions outside Vjv (again for a more precise definition 
see Section 2). The convergence result is the following: 


Theorem 2. Let Vjsj be as above and (i/’fl;)^g^rf be a DGFF with zero boundary conditions outside 
Vj^ with law Py^. Let the centering and scaling be as in (1). Then for a// z G R 


Ita p 

N—>+00 y fljy 



exp(—e ^). 


The core of the proof is an application of Slepian's Lemma and a re-run of the proof of 
Theorem 1. 

The structure of the article is as follows. In Section 2 we recall the main facts on the DGFF 
that will be used in Section 3 to prove the main theorem. 


^1 A| denotes the cardinality of A. 
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2. Preliminaries on the DGFF 

Let d > 3 and denote with || • || the £oo-norm on the lattice. Let rp = be a 

discrete Gaussian Free Field with zero boundary conditions outside A C . On the 

'—■ 

space n := R endowed with its product topology, its law Pa can be explicitly written 
as 

PA(di/^) = ^exp (JJ dQ{d%). 

\ ||a-^||=l JoceA aeZ'^\A 

In other words ipcc = 0 Pa"^- s. if a G \A, and {ipcc)aeA is ^ multivariate Gaussian ran¬ 
dom variable with mean zero and covariance {gA{‘>^> /l))^ ^ez‘^' gA is the Green's 

function of the discrete Laplacian problem with Dirichlet boundary conditions outside A. 
For a thorough review on the model the reader can refer for example to Sznitman (2012). It 
is known (Georgii, 1988, Ghapter 13) that the finite-volume measure ip admits an infinite- 
volume limit as A "I" Z'^ in the weak topology of probability measures. This field will be 
denoted as cp = It is a centered Gaussian field with covariance matrix g{Di, /3) 

for A, j6 G Z'^. With a slight abuse of notation, we write g{oi — (i) for g{0, a — j6) and 
also gA(^) = gAi^f ^)- It will be convenient for us to view g through its random walk 
representation: if Pa denotes the law of a simple random walk S started at a G Z'^, then 


g{a, p>) = Ea 


n>0 


In particular this gives g(0) < -|-oo for d > 3. 

A key fact for the Gaussian Free Field is its spatial Markov property, which will be used 
in the paper. The proof of the following Lemma can be found in Rodriguez and Sznitman 
(2013, Lemma 1.2). 


Lemma 3 (Markov property of the Gaussian Free Field). Let (D ^ K U := Z^ \R 

and define 

(Pcc = Ha, A G Z^ 

where ya is the (T{(pp, f G K)-measurable maip defined as 

^ Pa (Hk < +00, Shk = f) A G Z^^. (2) 

peK 

Here := inf {n > 0 : S„ G K} . Then, under P, (^a)^g^rf is indeyendent of (T{(pp, /3 G K) 
and distributed as under P(j. 


b 


A (E 


B means that A is a finite subset of B. 
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As an immediate consequence of the Lemma (see Rodriguez and Sznitman (2013, Re¬ 
mark 1.3)) 

p ^ = Pti {{% + e •) P-a.s. 

where is given in (2), P^i does not act on and has the law Pli- 


2.0.1. Law of large numbers of the recentered maximum. Although this can be obtained di¬ 
rectly by Theorem 1, we think it is interesting to insert an independent proof of the be¬ 
havior of the maximum of the DGFR 


Proposition 4 (LLN for the maximum). Let := 
following limit holds: 

E [max^gv^ cpc] 
N^+co -y/2 log N 


[0,n-lYnXf N 


= g{0)- 


n‘^ > 0. The 


Proof Observe first that g(0) > 1 (Lawler, 1991, Exercise 1.5.7). The upper bound follows 
from Talagrand (2003, Prop. 1.1. 3) with t := g(0) and M := N. As for the lower bound, 
we will use Sudakov-Fernique inequality (Adler and Taylor, 2007, Theorem 2. 2. 3). We 


first need a lower bound for d(x, f) 


2 

— cp^) : we will apply here the bound 





all > d 


( 3 ) 


whose proof is provided in Sznitman (2011). The key to obtain the result is to use a 

diluted version of the DGFF as follows. Gonsider := ^ where k := [log n\ G 

{1, 2, ...}. Note the fact that the expected maximum on is lower bounded by that on 

the diluted lattice Now for a, f ^ T := and k > d 


(3) 


« = Pm-2gw-f.) > 72 

> 72, j{0) - 


\ 


x(0) 


1 


I'/d 


i-2 


[lognj 


=: v(n, d)>0 


[Aif 


for n large enough. Notice also that limjv^.+oo r(iI/ d) = ^/2g(^. Hence by an application 
of the Sudakov-Fernique inequality 


E [max^gv^ cp^] 
V^ogN 


> v{n, d) 
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We obtain log |r| = rflog[^J (1 + o (1)) = rflog 
1 + o (1) and 


[logMj 


(1 + 0 (1))'^. It follows that — 


logN 


lim > ,/2s{0). 

N —>+00 Y^log A 1 


□ 


3. Proof of the main result 

The proof of the main result is an application of the Stein-Chen method. To keep the 
article self contained we recall the result from Arratia et al. (1989). 

3.1. Poisson approximation for extremes of random variables. The main tool we will 
use relies on a two-moment condition to determine the convergence of the number of 
exceedances for a sequence of random variables. Let {Xa)a^A be a sequence of (possibly 
dependent) Bernoulli random variables of parameter Pq.. Let W := ^ 

E [W]. Now for each a we assume the existence of a subsef Boi ^ A which we consider 
a "neighborhood" of dependence for the variable X^, such that Xa is nearly independent 
from Xp if j6 G A \ 5^. Set 

:= 5^ VocVfir 

b2-.= Yl E E[X.X^], 

aeA oi^^eBct 

&3:=5^E[|E[X«-p,|7^,]|] 

where 

'Hoc '■= cr (X^ : /3 G A \ Ba) • 

Theorem 5 (Theorem 1, Arratia et al. (1989)). Let Zhea Poisson random variable with E [Z] = 
A and let djv denote the total variation distance between probability measures. Then 

drvmw), CiZ)) <2{bi + b2 + b3) 
cttid 

P(W = 0) — e~^ < min|l, (^i + b 2 + b^). 

Let now A (s with N := \ A\, U'si(z) := a^^z + b]\j, and define for all a G A 

Mp^)- 

A standard tool to determine the asymptotic of p is Mills ratio: 




V^t 


< P (A/'(0, 1) > f) < 


e-^ 

Alnt' 


t > 0. 


(4) 


‘^/(N) = o (1) means limjv^+oo/(N) = 0. 
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This yields ^ N~^ exp(— z)^^. Since p^ is independent of a, we suppress the subscript 
a throughout. We furthermore introduce W := ^ [^] ^ Of 

course W is closely related to the maximum since {max^^gy^ cp^ < Ujv(z)} = {W = 0}. We 
will now fix z G R and A := e~^. We are now ready to prove our main result. 

Proof. Our main idea is to apply Theorem 5. The proof will first show that the limit is 
Gumbel, and in the second part we will prove uniform convergence. To this scope we 
define, for a fixed but small e > 0, 

(logN)2+2^) n A 

where B(a., L) denotes the ball of center a of radius L in the £oo-distance. We draw below 
examples of such neighborhoods when a G 9 A := | 7 GA: 3/3g \A, ||/3 — 7 || = l| 

and a. G int(A) = A \ 9A. 


A 

A\B^ 



Ba 



Z‘‘ 



(a) Bfl. when a G mt(A). 


(b) Ba when a G dA. 


Figure 1. Examples of Bft. 


~ y means that lim]v^+oo/(M)/y(N) = 1. 
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Convergence. The method is based on the estimate of three terms (cf. Subsec. 3.1). 

(i) Recall bi = Yl,oieA V^- Using Mills ratio we have 

I V27TOjv(z) 

= = o(l). (5) 

(ii) Recall h 2 = Yl,ix.eA ^ • First we need to estimate the joint probability 

P [cpa > Uf^{z),(p^ > Un{z)) . 

Denote the covariance matrix 


E2 = 


y(0) g{oc-^)' 

g{(x - g(0) 


Note that, for w G IR^, one has 


Iv’—1 

W E2 W = 


g ( 0 ) 2 _ g (^_^)2 b(0) (“-? + ”> 1 ) - 2^(“ - / i )“’ l “’ 2 ) . 


Using 1 := (1,1)^ we denote by 


A ; := Un(z) l^E 


I ^v —1 


UNiz){g{0) - g{ix - ^)) 


Wn(z) 


i = 1, 2. 


'2 1- ^(0) 

Exploiting an easy upper bound on bi-variate Gaussian tails (see Savage (1962)) we have 


P(<p. > un(z),<Pp > un(z)) < 


( 6 ) 


Note that using the explicit formula for the determinant one can bound the first factor 
easily by 


1 


1 


27r|detE2|U2AiA2 


< 


1 + 


§(o) ) 


1 _ gii^-b) \ 

g{Q) ) 


1 / 2 ' 


Now using M^f(z)2 = b^ + 2g{0)z + g{0)^z^/b^ and the bound of bfj (Hall, 1982, Equation 

3) 

g(0)(2 log N — log log N — log 47r) < bj^ < 2g(0) log N 
we can now upper bound the exponential term by. 


exp 


Mn(z)2 ^ 
2 


2 , \ ___ 2g(0)z 

^E2 ^1 1 < N g(o)+s(«-^) e s(o)+g(‘^-^)^ ^ '' 
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Also note that for X 7 ^ 0/S'(||'^||)/S'(0) < S'(^i)/S'(0) = 1 — KwhereK := Pq 
( 0,1) and Hq = inf {n ^ 1 : Sn — 0 }. Hence we have that 


= +00 G 


^(0) 


^(0)+^(n-^) 

We obtain thus 


> 


and 




2-k g{0)+g{‘^-^) 

(Z-'cP'C,- 


< 1-K. 


P [cpoc > UM{z),(pfi > um{z)) < - ^—N ( 2 -k) jnax (e 2^1{z<o},e ~^{z>0}) • 


We get finally for some constants c, c' > 0 depending only on d and k 

b 2 < cN(logN)^(^+^^) ~ (^ max (e-^" l{,<o},e- 2 "/( 2 -K) ^ 

< c'N (2-’^) (logN)^^^"^^^^^ max 1 

Since k/(2 — k) >0 we have that b 2 = o (1). 


{z>0} 


4z>0} 


(7) 


(hi) Recall b^, = ^ [l^ “ P« I l]- convenient to introduce also 

another cr-algebra which strictly contains Hoc = cr (X^ ; /3 G A \ Ba), that is 

'^'ix = ^ {(Pb- jS G A\Ba) . 

Using the tower property of the conditional expectation and Jensen's inequality 

E[|E[X«-p|it^^]|] <E[|E[X«-p|K]|]- 

At this point we recognize, thanks to Corollary 3, that 

E [Xa I n'^] = Pzrf\(A\B,)(V^'^ + > i^n(z)) P -a. s. 

where is a Gaussian Free Field with zero boundary conditions outside A\Boc. 

In addition, observe that A g(0) (Lawler, 1991, Section 1.5). We will write more 

compactly Ua := \ (A \ B^). 

We will make use of the fact that p^ is a centered Gaussian, and apply the same es¬ 
timates of Popov and Rath (2015): first observe using strong Markov property we have 
|6 G A \ B«, 

g{x,^) = ^ < +^r = 7 ) gil^^)- ( 8 ) 

7€A\B« 

We can plug this in to obtain 

Var [floe] = Y1 {^a\b, < gioc, < sup g{oc, 

beA\Bci b^A\Bci 

c 


< 


(log N) 2 (l+e)(d- 2 ) 


( 9 ) 
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by the standard estimates for the Green's function 

Cdik - jSfjS) < q||a - (10) 

for some 0 < Q < Q < +00 independent of a and fi (Lawler, 1991, Theorem 1. 5. 4). 
Using the estimate 

P(|A/'(0,1)| > fl) < fl > 0 (11) 

we get that there exists a constant C > 0 such that 

p {M > (rn(z))“^“") < Cexp (-(logN)(2‘^-5)(i+")) . (12) 

Note that this quantity goes to zero since d > 3. Hence 


E 
+ E 


+ > Um{z)) - p 

PuA^a + Pa > U^iz)) -p 


= E 

1 


+ Pa > Un{z)) - p 




[\}ia\>{UNiz)) ^ ‘"I 


= : Ti + T2. 


By (12) and the fact that d >3, we notice that NT 2 = o (1). Therefore if is sufficient to treat 
the term Ti. By conditioning on whether p is larger or smaller than Pn^ (ipa. + Pa > Wn(z) ) 
we can split the event in Ti into the following two terms. 


E + Un{z)) - p) l{,<p„.(»,+„.>„„(.))}_ 

+E [(p - PuTf. + P. > Mz))) 1{,>P„.(^+„.>.„W)} 

=: Ti^i + Ti^2- 


(13) 


We will now deal with Ti 2 - The first one can be treated with a similar calculation. Using 
Mill's ratio and fact that xpoc has variance gu^ we get that. 


P - PuA% + Pa > Un{z)) 

^ _ L_f VguM) 

~ I yWN(z) — fla 


(upi(z)-l‘aP 

VguM)^ 

^/iniuNiz) - Poc) 


We have on the event 


\Pa\ < {un{z)) 


— 1—e 


that the above is bounded by 


(14) 


\/2nuj^(z) 


/ 

l-(l + o(l)) 

V 


[I g ( 0 ) ^ » n ( 9 ^ I » Af ( z ) 

VguMuN{z)e\ 

^/g(U)UN{z){l - Un(z)-^-^) 


»n(^) ^ \ 

(“) 


/ 

( 15 ) 
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Since the bound is non random, by bounding the indicator functions by 1, 
E 

Now 


(p - Pu.(^« + > «n(z)) j .<„„«)} 


< (15). 


( 12 ) 


&3<5^(ri + T2) < 5 ^Ti + o(1) = 5^Ti4 + 5^Ti,2 + o(1). (16) 

aeA aeA aeA aeA 


Then 


V^e 2,(0) 


/ 


Ti,2 — -^ 


\^2nuN{ 


/ 


l-(l + o(l)) 


V 


VSuMun{z) 


1 _g(o) ipnA 

2?(0) ^ ^ 


V 


v/^mn(z)(1 + o(1)) 


/ 






Observe that 1-< 0 since g{0) > gUa(^)- observe further that (and we will 

prove it in a moment) 

Claim 6. sup^g^ (l - Uj^{zY = o (1). 

Therefore Ti 2 = o (1) uniformly in a. This yields that 

^ Tp 2 < (1) = o (1). (17) 

V27TOn(z) 

Analogously, Yl,ix.eA ~ ‘^ (!)• (1^) ^ (16), one obtains b^, = o (1). We now 

only need to show Claim 6. By the Markov property we know 

gUx{‘^) = ^( 0 ) - {Ha\Bx < +^' = 7 ) 

7€A\B„ 

This shows that 

0 < g(°) _ 1 < "“P7e.4\BA(0'-‘«) 


^LI,(a) gUx[‘>^) 

Note thatg( 7 ,a) V Cd(logN)-2(^-2)(i+e). Also,gtj„(a) = 1{S„=«} 

hence we have 

,(0) 1 ^ .n^,.Ml-2(d-2)(l+e) 


0 < -1 < c(logN)- 


gUx[‘^ 


> 1 and 
(18) 


from which it follows that 
^(0) 


fl- UN{zf < c(logN)-2(^-2)(i+^)(logN + z + o(l)) = 0 ( 1 ). (19) 

V gUx{‘^)J 


Therefore the claim follows and we have shown pointwise convergence. 


□ 
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3.2. DGFF with boundary conditions: proof of Theorem 2. The idea of the proof is to 
exploit the convergence we have obtained in the previous section. We will show a lower 
bound through a comparison with i. i. d. variables, and an upper bound by considering 
the maximum restricted to the bulk of Wf, concluding by means of a convergence-of- 
types result. We abbreviate gN{', •) ■= ^ > 0 define (recall that Wf = 

[0, n - 1]^^ n with N = 

< := |a G % : ||a - 7ll > 7 G \Vn } . 

We begin with the easier lower bound. 


Proof of Theorem 2: lower bound. We will need a lower and an upper bound on the limiting 
distribution of the maximum. Let us start with the former. We use the shortcut Pjsi := Py^. 
First we note that since the covariance of (xpoi) is non-negative, we can apply Slepian's 
lemma for the lower bound. Let {Za)aeVj.j be independent mean zero Gaussian variables 
with variance then by Slepian's lemma it follows that 


Pn 


( max Za < un(z) 
\oceVN 


< Pn 


( maxipc < Un(z) I / 


where Ufq{z) = a^z + &n as before. Then we want to analyze P{m.aXi^(^j^Za, < u^{z)). 
First fix z G R. Take N large enough such that —g{0)bfi < z (this is possible as > -|-oo). 
Now note that 


Pn 



(4) 
> 




e 

\/2nuM(z) 



> 


/ 

1 

V 


»Nb)^ 
e 2?(o) 

^/2nu]s|(z) 


\ 

s/m 

/ 


N 


The last term converges to exp(— e as N m + 00 . This shows that for any fixed z G R, 


liminf P^ 

N—>■+00 


( maxxpa < U]^{z) ] > exp(—e 

KocEVn 


□ 


In order to prove the upper bound of Theorem 2, we shall need a Lemma which will 
allow us to derive the convergence of the maximum in V/v from that of the maximum in 

^N- 


Lemma 7. Let N > l,T^be a distribution function, and = (1 — 2d)‘^N. Let and b^q be 
as in (1). 7/lim]v^+oo fN(«mNZ + bm^) = exp(- e“^), then 


lim Ffq{afqz + b]si) 

N—>+00 


exp (^_e-z+'^iog(i- 2 J)j _ 
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Proof. The proof follows from a convergence-of-types theorem (see Resnick (1987, Propo¬ 
sition 0.2)) if we can show that 


On 


1 and 


Un 


dlog(l — 2^). 


It is easy to see that 


^m]v 

On 


^ I dlog(l-2^) \ 


1/2 


logN 

To show the second asymptotics note that 

dlog(l — 2S) 


1 . 


/2g{0) log mw - y 2^(0) log N = 
Also observe that as N 


+ 0 


m 


21 ogN ' “ V(logN)^ 

-|-oo one gets 

loglog(47rN) loglog(47rmN) 
2y/2To^ 2 a/ 2 logmjv 




- log 1 + 


(ilog(l — 2S) 


2V'2Tdpi 

So using the above equation and (21) we get that 

t’ffXjAf ~ t’N _ ~ b]\j 


fljV 


m 


logN 

\/2gi0) logN(l + o(l)) 


( 20 ) 


2^(0) logN. (21) 


+ 0 ( 1 ) 


(ilog(l — 2S). 


□ 


We have now the tools to finish with the upper bound. 

Proof of Theorem 2: upper bound. First fix z G R and ^ > 0, set mw := \V^\ = (1 - 2^)‘^N. 
For the upper bound, we again use Lemma 3 and the fact that, for a G V^, one has the 

equality (pa = + Pa^'^ under the infinite volume measure P, where = E [^^1 . 

Hence if we fix e > 0, and condition on the event that 

max ^\ < eum^ 

(where flmjv is defined according to (1)) we have 


P]V ( max ^a < 

\aeVN 


u 




z) I < P max (pa < Ufn^{z -E e) j -E P max \ pa\ > 






niN 


( 22 ) 


First we show that the second term goes to zero. Observe that pa is a centered Gaussian 


with variance 

maxVar [p^] < sup g(/3, 7 ) = O 

pev^,jedVj^ ^ 
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Let be a collection of i.i.d. Gaussians with mean zero and E = E for 

all A. By Slepian and Talagrand (2003, Prop. 1.1. 3) we have 


2E 


P max \}ia\ > < 








+ 2 P max Oft. < 0 


. o^evi 


N 


< 


i^max^g^, Var log |Vj| 




+ 0 ( 1 ). 


Since fljv grows like (a/ 2logN) ^ as N —> +oo, we can conclude that, for every e > 0, 


lim P 

N—> + 00 


max |?ift:| > eumj^ 


Using Theorem 1 we have that 


0 . 


lim P max cpa. < Umj^iz + e) ) = exp ( — e 

N->+oo yaev^ y 

and, being the limit continuous, we let e —^ 0 obtaining from (22) 


limsupP^/ maxxpa < Umj^{z) < exp(—e ^). (23) 

N—>+00 J 

Now using an easy comparison with independent random variables just as in the proof 
of the lower bound of Theorem 2 above it follows that (23) is in fact an equality. By 
Lemma 7 one can conclude that 


Pn 



exp (_e-^+'^iog(i- 2 + 


and thus letting ^ 0, the upper bound follows. 


□ 
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